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Abstract 

This paper describes the procedures for and 
results of aeroservothermoelastic studies. The objectives 
of these studies were to develop the necessary procedures 
for performing an aeroelastic analysis of an 
aerodynamically heated vehicle and to analyze a 
configuration in the classical "cold" state and in a "hot" 
state. Major tasks include the development of the 
structural and aerodynamic models, open loop analyses, 
design of active control laws for improving dynamic 
responses and analyses of the closed loop vehicles. The 
analyses performed focused on flutter speed calculations, 
short period eigenvalue trends and statistical analyses of 
the vehicle response to controls and turbulence. 
Improving the ride quality of the vehicle and raising the 
flutter boundary of the aerodynamically-heated vehicle up 
to that of the cold vehicle were the objectives of the 
control law design investigations. 
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Introduction 

The proposed flight envelope of the National 
AeroSpace Plane (NASP) includes Mach numbers from 
low subsonic to hypersonic speeds. As a result, the flight 
vehicle will experience a wide range of dynamic pressures 
and aerodynamic heat loads which may degrade the 
aeroelastic response characteristics of the vehicle. The 
single-stage-to-orbit mission of the NASP requires a large 
vehicle with a high fuel mass fraction. These conditions 
generally result in flexible configurations having low 
structural frequencies. Thus, both the flight environment 
and the required vehicle characteristics may cause 
structural designs to exhibit aeroelastic instabilities. 

The work described in this paper is part of a 
research effort at the NASA Langley Research Center to 
develop analytical methods for predicting the aeroelastic 
characteristics of hypersonic aircraft with active control 
systems. Major tasks associated with this activity 
include: 1) construction of a generic hypersonic model for 
use in performing these analytical studies, 2) development 
and analysis of aeroelastic models which incorporate the 
effects of aerodynamic heating, and 3) application of 
active controls to improve aeroelastic response. Flutter 
suppression and ride quality augmentation systems were 



designed and evaluated for the aircraft in its classical cold 
state and for the aerodynamically heated state. 


Modeling and Equations of Motion 

Aeross iw i temire lastic Modeling 

Aeroservothermoelastic modeling, by definition, 
incorporates aerodynamics, controls, temperature- 
dependent material properties and the structural definition 
of the configuration. The flow chart in Figure 1 presents 
the analysis and control system design approach taken. 
The three major aspects of the model development are 
described in the boxes on the left. These include 
determination of thermal loads, construction of a finite 
element model and calculation of unsteady aerodynamic 
forces. These data are combined to form aeroelastic 
equations of motion for the aerodynamically-heated 
vehicle. Open loop analyses are performed on these 
models which also serve as a testbed for control law 
design. After control laws have been formulated, closed 
loop analyses and simulation may be conducted. 

A more detailed accounting of the procedure for 
the model development task [1] is given in Figure 2. The 
individual processes are indicated by bracketed numbers. 
The initial step is to perform a trim analysis (1) of the 
vehicle in a rigid and unheated state to find the 
approximate flight attitude and control surface settings. A 
steady aerodynamic analysis (2) is then performed to 
determine the external pressures, the aerodynamic heating 
rates and the radiation equilibrium wall temperatures. 
These radiation equilibrium wall temperatures are used to 
approximate the temperature loading on the elastic finite 
element model. A thermal structural analysis (3) is 
performed to determine the temperatures at each of the 
node locations. New vehicle shape and stiffnesses are 
determined by performing a static structural elastic 
analysis (4&5). The nodal temperature data and the 
temperature-dependent material properties are incorporated 
with the geometric and element definitions to yield the 
finite element data sets which contain structural properties 
for each element based on the local temperature. The 
temperature-induced internal stresses on the vehicle are 
also computed using the nodal temperature data as a 
thermal load case. These stresses are incorporated into the 
stiffness matrix, which has already been changed to reflect 
the material properties changes associated with heating the 
vehicle. A classical vibration analysis (6) is then 
performed using the effective stiffness matrix. The 
resultant mode shapes and natural frequencies are input 
into an unsteady aerodynamic analysis (7&8) to determine 
generalized aerodynamic forces.(GAFs), required for 
stability analyses, control law design tasks, etc. 

Geometric and Finite Element Model Descriptions 

Isotropic Titanium- Aluminide was used as the 
material model for the generic configuration. The 
temperature dependence of material properties, (the 


Young's modulus, the shear modulus and the thermal 
expansion coefficient), are included in the modeling. The 
geometric and inertial properties of the configuration are 
detailed in Figure 3 and Table 1. The finite element 
model, shown in Figure 4, consists of five components: 
the fuselage, wingbox, aileron, flap and wing tip-mounted 
fin. The fuselage has an elliptical cross-section with the 
width being twice the height. Ring and skin 
construction, using very stiff rings, prevents shell or 
"breathing" modes of the fuselage. The low-aspect ratio 
wing is constructed with spars, ribs and skins, and is fully 
attached to the fuselage. The wing has a 70-degree leading 
edge sweep and a 3-percent circular arc airfoil. The 
control surfaces consist of an all-movable vertical fin, 
flaps, and ailerons. 

vibration Mods Shapes 

Two sets of structural vibration data (mode 
shapes and frequencies) were generated: one which 
contained the aerodynamic heating effects and one which 
did not. All heating effects were based on the Mach 4 
flight condition. Both structural models, cold and hot, 
were analyzed at Mach 2 and Mach 4. Thus, the hot 
Mach 2 flight condition was assumed to be representative 
of the deceleration phase of the mission profile. The 
Mach 4 cold model was assumed to be in an acceleration 
phase such that the vehicle has not yet been affected by 
the aerodynamic heating. Since there were no visible 
differences between the "hot" mode shapes and the "cold" 
mode shapes, only one set of vibration mode shapes is 
shown in Figure 5. Note that the first and third elastic 
modes include significant fuselage motion. The 
frequencies, hot and cold, for each of the first six modes 
are listed in Table 2. As previously mentioned, the 
mission requirements cause the frequencies for this class 
of vehicle to be very low. Without aerothermal effects, 
the frequencies lie between 3 and 1 1 Hertz. When the 
effects of heating were incorporated, the frequencies drop 
by 13 to 20 percent, indicating the destiffening effect of 
the aerodynamic heat loads. 

Aerodmmigs 

Supersonic generalized aerodynamic forces were 
generated using van Dyke's second order Piston Theory 
[2], including the effects of airfoil thickness and non-rigid 
chords. The validity of the calculated GAFs was 
substantiated by computing stability derivatives and short 
period behavior and comparing the calculations with 
steady-state results [3]. 

Piston Theory. - The premise behind Piston 
Theory aerodynamics is that the relationship between the 
local pressure and the normal free-stream velocity at a 
point is analogous to the relationship between the 
pressure in a one-dimensional piston chamber and the 
velocity of the piston. At sufficiently high Mach 
numbers this analogy is valid, and local wave theory 
gives a good approximation of the unsteady aerodynamics. 
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The local, linearized pressure equation used in this Piston 
Theory Method is given as: 


when transformed into the Laplace domain the equations 
of motion can be written: 


Ap(x,y,t) = -2pa 


l + G^Z(x,y) 
dx 


w i 



where 

c M 4 (y + 1)-4P 2 
2 (3 2 
and 

P=Vm 2 -i 


The variable "z" represents both the spatial dependence and 
the time- varying displacement of a point. This displace- 
ment is calculated at discrete points on the surface of the 
vehicle by representing the wing surface and projection of 
the fuselage onto the x-y plane with trapezoidal boxes. 
Using surface spline interpolation the displacement of the 
center of the boxes was computed from the mode shapes. 
Normal velocities over the wing and fuselage projection 
surfaces were computed by assuming simple harmonic 
motion of each mode. The upper case "Z" is the distance 
from the mean camber line to points on the skin, thus 
representing the thickness. The circular arc airfoil 
thickness characteristics of the wing were incorporated 
into the pressure equation, however, the thickness effect 
of the fuselage was not included in the aerodynamic 

modeling.* The point forces resulting from the Piston 
Theory pressures were also concentrated at the center of 
each box. The generalized aerodynamic force for each 
mode was generated by summing these point forces, 
weighted by the interpolated mode shapes, over the wing 
and fuselage surfaces. 

Equations of Motion 

Classical open loop aeroelastic equations of 
motion were formulated for the longitudinal degrees of 
freedom. These equations of motion, based on two rigid 
body and the first six vehicle vibration modes are derived 
using Lagrangian energy equations. They represent the 
summation of inertial, dissipative and internal restorative 
forces along with the reduced-frequency dependent 
aerodynamic forces which are induced by the rigid body 
and structural motion, control deflections and gusts. 
Small perturbation equations, referenced to level 
equilibrium flight, were constructed to represent a 
linearized system of equations expressed in terms of 
generalized coordinates. As described in Reference 4, 


* At the time that the analysis and control law design 
were performed, the Piston Theory computer code did not 
have the capability of including the thickness effects of 
the fuselage. The assumption that the contribution to the 
GAFs would cause only small changes in the dynamics 
was made and subsequently validated. After enhancing the 
Piston Theory code to incorporate the fuselage thickness 
effects, the flutter dynamic pressure increased by two 
percent. 


([M]s 2 +[D] S+ [K])5 + q[Q|j 5 5 J = [Q,]^- 

where £ are generalized displacements of the rigid and 
flexible modes. 

Aerodynamic Approximations 

Generalized Aerodynamic Forces (GAF's) 
predicted using classical unsteady aerodynamics codes 
result in matrix sets of data at discrete values of reduced 
frequency. In order to incorporate the data into the 
Laplace domain aeroelastic equations, they must be 
transformed into continuous functions of the Laplace 
variable s. To this end, an approximation is generally 
made, employing a curve fit to the tabular data [5]. By 
examining a typical variation in GAF with reduced 
frequency, as seen in Figure 6, it is observed that, in, the 
case of Piston Theory aerodynamics, the real parts of the 
GAFs are constant and the imaginary parts are linear with 
reduced frequency. This characteristic allows the Piston 
Theory aerodynamics to be represented exactly by a first 
order equation: 




= Aq+Aj-s 


where Ao is the constant value of the real part and Ai is 
the slope of the imaginary part. 

State Snace Equations 


Control system design and analysis are more 
readily performed using a state space representation of the 
equations of motion. The transformation of the equations 
into first order form results in 


x= Ax+Bu , 

where 
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Flu tter 


as described in detail in Reference 4. Note that the 
aerodynamics, when expressed by the Ao and Ai 
coefficients, appear in the linear time-invariant equations 
as damping and stiffness increments. The above 
equations describe the flexible airplane dynamics. They 
were augmented with equations representing two third- 
order actuators, which brought the total number of states 
to 22. The output equations consist of angle of attack, 
pitch rate of the center of gravity and the pilot station, and 
the normal acceleration at six fuselage and wing locations. 


Qgcn Lm Analyses 

To identify the open loop plant characteristics, 
trends in the rigid body and aeroelastic characteristics with 
Mach number, aerodynamic heating and structural 
flexibility were calculated [6J. 

Short Period Dynamics 

The short period mode dynamics are influenced 
by both the aircraft flexibility and by the aerodynamics. 
Figure 7 provides a root locus plot which shows the 
changes in short period behavior incurred due to the effects 
of structural flexibility, aerodynamic heating and Mach 
number. The six curves represent the trace of the short 
period eigenvalue in the complex plane as the altitude is 
varied (altitude decreases from right to left along each 
curve). 

In an aeroelastic system, the roots of any one 
mode are influenced by the other modes near it Because 
the structural frequencies for these configurations are low, 
in the neighborhood of the rigid body frequency, it is 
anticipated that they would exhibit a large degree of 
influence over the short period mode. This influence can 
be seen by examining the roots of the rigid vehicle versus 
the eigenvalues after the effects of flexibility have been 
included. The figure indicates that including flexibility 
tends to have a destabilizing effect, since the roots move 
closer to the imaginary axis, which is the neutral stability 
boundary. The effects due to the aerodynamic heating can 
be evaluated by comparing the hot and cold data for the 
same Mach number. At either Mach number, the 
destabilizing effect of the heating is seen as the roots for 
the hot data fall further to the right in the s-piane than 
those corresponding to the cold data. 

To determine the effects of the Mach number, the 
curves for the rigid data, the hot data and the cold data 
must be examined separately. It is seen that as the Mach 
number is increased, the short period frequency is 
increased and the damping is decreased. Thus, increasing 
Mach number also has a destabilizing effect on the short 
period dynamics. 

Comparing the curves in these three ways shows 
clearly that the Mach number has a much larger influence 
than either the flexibility or the heating for the condition 
analyzed. 


The altitude root locus plot of Figure 8 shows 
some of the open loop flutter characteristics. These 
results are for the heated model at Mach 2. The 
eigenvalues associated with the vibrational modes are 
plotted as altitude is lowered. Since this is a matched 
point analysis, both the density and the velocity are 
changed. The predicted flutter mechanism involves the 
coalescence of the third and fourth elastic modes at a 
dynamic pressure of 4800 psf and at a frequency of 6.7 
Hertz. Figure 9 shows the predicted hot and cold flutter 
boundaries for this vehicle with respect to Mach number. 
The dynamic pressure was normalized by a reference 
dynamic pressure, the flutter dynamic pressure at Mach 4 
without heating. The region below either curve represents 
the area of flutter-free flight for that configuration. As the 
Mach number is increased, the models require higher 
dynamic pressures be applied before they become 
unstable. Heating the model lowers the flutter boundary 
over the entire range of Mach numbers examined in. this 
analysis at Mach 4, the drop in flutter dynamic pressure 
due to aerodynamic heating was approximately 25 percent. 

Frequency Response 

The state space equations of motion were 
augmented with a second order Dryden gust model and a 
first order roll-off filter to analyze the gust responses of 
the system. Open loop transfer functions between the 
gust input and the normal accelerations on the fuselage 
indicated dominance of the first and third flexible modes 
for these outputs. These two modes were previously 
shown to exhibit substantial fuselage motion. The wing 
bending and wing torsion modes were nearly unobservable 
from the pilot station acceleration due to gust input as 
shown in the power spectral density plot of Figure 10. 
This figure also illustrates the shift in frequencies due to 
aerodynamic heating. Although not shown, open loop 
transfer functions between the control surfaces and the 
pilot station acceleration were also examined as indicators 
of controllability. 


Active Control Laws 

One of the objectives of this project was to 
assess the application of active controls technology to 
counter dynamic responses aggravated by various flight 
characteristics. As seen in the open loop analyses, 
aerodynamic heating lowers the flutter boundary 
substantially. Thus, the first control concept investigated 
was flutter suppression. The intent was to recover the 
flutter dynamic pressure lost due to heating; that is, to 
raise the flutter boundary back up to its original, unheated 
boundary. The second active control concept investigated 
was a ride quality augmentation system (RQAS). The 
length of the vehicle, the aft center-of-gravity location and 
the low structural stiffness imply that small perturbations 
will induce large oscillations at the cockpit The RQAS 
was designed to reduce the accelerations at the pilot 
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station Figure 1 1 shows the basic feedback scheme used 
for each of these independently-designed and evaluated 

control laws. 

Flutter Sup pir^ 00 System 

The objective of this task was to recover the 
flutter dynamic pressure reductions caused by aerodyn^ 
heating through the use of active controls 
was also to determine the sensitivity of the flutte 
suppression system (FSS) design to changes in flight 
conditions. 

Procedure . - A full-order state estimator 
fKalman Filter) was used for compensation in the 
feedback loop The controller was designed using standard 
Gaussim. (LOG) cmmol to * tst*. 
methods with Loop Transfer Recovery (LTR) [7] to 
improve stability robustness in the face of changes in 
flight dynamic pressure. Normal acceleration at the pilot 
station and a t a location very near the wing aiieronwere 
used as measurements for feedback to the compensator 
Both measurements were assumed to be ^y. m FSS 
control law was designed to minimize the total system 
energy by weighting the sum of the structural strain and 
S energies aSd ^ comm^d^ control 
deflections (a measure of control energy). Tne g 
SI design was Mach 2 jus. beiow the Hauer 
dynamic pressure using the hot vehicle dynamics. 

The FSS control law is described by the 
equations 

u = Fxf 

Xf = Axf + Bu + H(y-Cxf) 

where the gain matrices F and H come from the solutions 
of the linear regulator and Kalman Filter pro cms. 
dynamics for the estimated states xf are appended to th 
site-space dynamics of the open loop aeroservothermo- 
elastic aircraftmodel for closed loop performance analysis. 


Results . - The performance of the FSS control 
law in recovering the cold vehicle flutter dynamic pressure 
was evaluated by determination of the eigenvalues of th 
closed loop hot vehicle equations of motion at S 

Syndic pressures. Tbe dynamic pressure was mcreared 
bv reduction in flight altitude using the match point 

densftyfortoth Mach 2 and 4. Tbe dynamic pressme was 
increased until an unstable eigenvalue (indicating flutter) 

Zs obtained or until the altitude reached sea M. The 
results are shown in Figure 12 as the ratio of hotwe^d 
flutter dynamic pressure from Mach 2 to 4 for both the 
open and closed loop cases. The open loop data 
previously shown in Figure 9, is along the loww edge of 
die box with a ratio less than 1.0, ‘ndicaung die loss^n 
flutter dynamic pressure due to aerothermal heating^ The 
closed loop results are along the upper edge of the box, 
wSLe tat representing the increase in 
pressure due to the FSS control law. Note that the FSS 


increased the flutter dynamic pressure of the hot vehicle 
well beyond that of the cold vehicle, with the Mach 2 
result limited by meaningful altitude rather than actual 
flutter Note, too, that the results for both Mach 2 and 4 
“Tor Z same control law indicating the relanve 
robustness of the control law to variations in Mach 
number. This implies that a FSS system may be 
implementable on an actual vehicle without extensive 

gain scheduling. 

pH* Qualit y /Augmentation Svstem 

A ride quality augmentation system (RQAS) was 
designed to assess the potential for reducing the 
turbulence-induced normal accelerations fe t y e p 
due to excessive fuselage flexibility. Additional^, the 
influence of the aerothermal loads on the closed loop 
system were to be evaluated. 

pnas Philosophy . - The vehicle vibration mode 
shapes indicate that the fuselage bending modes pccur at 
K^Sea and iuHueuce the aircraft dynm*« »» 
large Stem. Open loop analyses have shown mat the 
normal acceleration response to turbulenceat the pi 
station is dominated by the fuselage tnmles "ap 
ailerons are located on the trailing edges of the wings, 
nearly the entire length of the vehicle away from the point 
nf interest. Without the addition of a canard, fuselag 
flexibility prohibits the collocation of the feedback sensor 
and control surface. The accelerometer measurement at 
the pilot station is used directly as the feedback signal, 
despite complications due to noncollocauon. 

n«i«m Model. - The ride quality augmentation 
system ^ designed at Mach 4 at 30000 feet, 
corresponding to a flight velocity of 3980 feet per second. 
Basedon statistical responses and transfer function data, 
Seated model was more affected by turbulence than the 
cold model; the indication being that more control ^ ) ^ 
ts required to compensate the performance of the hot 
vehide. Thus, the RQAS was developed based on the 
heated configuration and tested for robustness in the fac 
ifdre aerodynamic heat loads by applying the 
compensator to the unheated configuration. 


r* nntml I -aw Design - - The RQAS control law 
was constructed using the same sensor and noise 
intensities applied for the FSS; the controller was 
designed via an eigenvalue assignment t^hmque ffl. 
transfer function of primary interest m die RQAS design 
was the normal acceleration response at the pilot station 
dS te gust disturbances. The flap was used as the 
feedback control surface, after controllability tests 
indicated that it had several times more influence on this 
rSSSe than the aileron. In placing the eigenvalues 
much consideration was given to the fact that the natural 
freaSeS Se subject to relatively large variations 
brought about by the aerodynamic heating. Due to i these 
knoiJn "uncertainties," it was desirable to expendve^ 
Stie control effort altering the frequencies 
eigenvalue locations were chosen such that they pre 
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shifts along horizontal lines in the complex plane Th* 

WCre to produce a critically daSSd shw 
penod mode and to increase the damping of thefost and 
second bending modes of the fuselage^ ?i£r o^ 

tan, , ^ CS ?* t S ~ ' The Performance of the RQAS control 
law was evaluated by determining the reduction in the 
peak acceloations felt at the pilot station and also in the 
reduction of the root mean square (rms) resDonsTfV^r^ 
surface limitations were set at 10 degreesof deflection 
with maximum rates of 30 degrees per second Thee? 
limits were evaluated for a gust velocity of 41 68 feet per 

? .f 1,ght condltion (30,000 feet altitude! 
Additionally, the control law must not suppress the 

SS" lie 'm 0 * My " xwm, J ^ase il 

weiToSedV, "|fP onses “ vehicle localions 

criteria Th» M M mdlcator o{ ^ violation of this 
entena. The rms responses for the heated and unheated 

as operating as shown in Figure 13 The closed ir^« 

set,,*: k° 

a 53.8-percent^E' faSS lSSg 

for the closed loop heated vehicle and T nsSS 
reduction for the cold configuration. PCrCent 

Conchidinp p fmnr | n 

fh __ , Thls papcr has described the aeroservo 
modeling > analysis, and control law design 
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Arrows Indicate Increasing dynamic pressure 


MOOE 

DESCRIPTION 

FREQUENCY (HZ) 

% REDUCTION 



COLD 

HOT 

/COLD-HOT) 
' COLD ' 

1 

fuselage 

bending 

3.01 

2.43 

10.3 

2 

wing 

bending 

4.02 

3.46 

13.4 

3 

second 

fuselage 

bending 

7.06 

5.67 

10.7 

4 

wing 

torsion 

7.70 

6.56 

14.6 

5 

ail-movable 
fin mod# 

0.47 

7.63 

104 

6 

highly coupled 
mode 

10.06 

8.64 

10.3 


Table 2. Vibration Mode Descriptions and Natural 
Frequencies. 
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Figure 8. Altitude Root Locus for Mach 2 Heated 
Configuration. 
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Figure 6. Lift Due to Pitch Displacement Genralized 
Aerodynamic Force. 


Figure 9. Flutter Boundaries. 
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Figure 7. Effects of Mach Number, Heating, and 
Flexibility on the Short Period Mode. 
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Figure 10. Open Loop Normal Acceleration Response 
at the Pilot Station Due to a 1 ft/sec. Gust 








Commands 



Figure 1 1 . Schematic of Feedback Control System for 
FSS andRQAS. 



Figure 12. Flutter Suppression Results. 
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Figure 13. Power Spectral Densities: Normal 

Acceleration at the Pilot Station Due to 
a 1 fL/sec. Gust 
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loop vehicles. The analyses performed focused on flutter speed calculations, short period eigenvalue 
trends and statistical analyses of the vehicle response to controls and turbulence. Improving the ride 
quality of the vehicle and raising the flutter boundary of the aenxlynamically-heated vehicle up to that of 
the cold vehicle were the objectives of the control law design investigations. 
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